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Wave effects can be important for the gravitational lensing of gravitational waves. In such a 
case, wave optics must be used in stead of geometric optics. We consider a plane wave entering 
a lens object and solve numerically the wave equation for three lens models: the uniform density 
sphere, the singular isothermal sphere, and the Hernquist model. By comparing our numerical 
solutions with the analytical solutions under the thin lens approximation, we evaluate the error of 
this approximation. The results show that the relative error of the thin lens approximation is small 
if the geometrical thickness of the lens is much smaller than the distance between the lens and the 
observer. 
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I. INTRODUCTION 

Gravitational waves from coalescing compact binaries 
composed of neutron stars or black holes are the most 
promising targets for ground-based as well as space-based 
detectors. By applying the matched filtering technique 
which uses our theoretical predictions of wave forms ob- 
tained by the post-Newtonian computations, we can ex- 
tract the binary parameters such as the masses of each 
compact object, distance to the source, spatial positions 
of the source and so on [| . 

One possibility which alters the predicted wave forms 
calculated with high precision is the gravitational lensing 
of gravitational waves. If a massive object lies suitably 
between the source and the observer, gravitational lens- 
ing of gravitational waves occurs. One important point 
is that since the wavelength of gravitational waves we are 
interested in is much larger than that of light, a situation 
where the geometrical optics approximation breaks down 
can be realized in some cases. As is discussed by many 
authors 0, 0, 0, IE 1 if the wavelength is larger than 
the Schwarzschild radius of the lens object, the diffrac- 
tion effect becomes important and the magnification ap- 
proaches to unity. Therefore we must use wave optics 
rather than the geometric optics for 

M L < 10 8 M Q , (1) 

where Ml is a mass of the lens and / is the frequency 
of the gravitational waves. This frequency (mHz) is the 
case for the planned detector LISA J?}. 

Further because the gravitational waves from a com- 
pact binary are coherent, interference between lensed 
waves is important. Note that this situation is not in 
general realized in the case of gravitational lensing of 
electro-magnetic wave such as visible light. Since light 
is emitted from microscopic region (usually atomic size) 
which is much smaller than the size of the source, each 
phase of the electro-magnetic wave emitted from differ- 
ent points has no correlation and thus interference effect 



vanishes. 

If we assume the coalescence of SMBHs of mass 10 4 ~ 
10 7 Af Q as the source of the gravitational waves, it can 
be detected even if the sources are located at the cos- 
mological distance (z > 5). Event rate of SMBH-SMBH 
merger for LISA is estimated as 0.1 ~ 10 2 event/yr Q 
and lensing probability becomes several percent. Hence, 
some lensing events per year will be detected by LISA. 

Motivated by the fact that wave effects can be detected 
for the gravitational lensing of gravitational waves, there 
are now growing interests in the wave optics in gravita- 
tional lensing gESIllEIHQEm- However 
we have to solve wave equation which is generally partial 
differential equation between the source and the observer 
in order to obtain the lensed wave form at the observer. 
Except for a few special cases, exact solutions of wave 
equation are not known at present. Several authors have 
used thin lens approximation which reduces wave equa- 
tion to double integral for single lens object (for multi- 
lens objects, integration becomes multi-integral, see |l4^ 
and thus makes the problem easier In geometric 

optics the trajectory of light ray is obtained by solving 
geodesic equation, and it is known that thin lens approx- 
imation is valid However there has been no studies 
or comments about the validity of the thin lens approxi- 
mation in the framework of wave optics. 

In this paper, we develop a formulation to solve the 
wave equation for a spherically symmetric lens, where a 
partial differential equation reduces to a set of ordinary 
differential equations. We also solve those equations for 
simple lens models: the uniform density sphere, the sin- 
gular isothermal sphere and the Hernquist model and 
evaluate the error of the thin lens approximation. 

This paper is organized as follows. In section II, we 
briefly review wave optics in gravitational lensing under 
the thin lens approximation. In section III, we develop 
a formalulation to solve the scattering problem of grav- 
itational waves by a lens. In section IV, we present our 
numerical results and discuss the validity of the thin lens 
approximation. Section V is a summary. We use unit of 
c= 1. 
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density £(s), where s — (x,y). The two dimensional 
potential ?/>(s) in z = plane is defined as 



(4) 



ip(s) = 2 / dzU(s,z). 



Here, ip is also obtained from the surface density by using 
V^(«] = 8ttE(5). 

It is useful to define the amplification factor F (which 
is called the transmission factor in Ref.01) as F — 4>l4>oi 
where <f> is the gravitationally lensed waveform obtained 
by solving Eq.Q and </>° is the incident wave. The am- 
plification factor at the observer under the thin lens ap- 
proximation is given by |l7j | 



F t hin(r,9) 



2nir 



incident wave 



(f se l ^\ r9 - s \ -^V-(s). (5) 
Here, F t ^ in is normalized so that F t ^ in = 1 for ip = 0. 



FIG. 1: Lensing configuration. The lens is the origin of coordi- 
nate axes, while the observer position is (r, 6) with 9 < 1. The 
incident wave is a plane wave propagating in the z-direction. 

II. GRAVITATIONALLY LENSED WAVEFORM 
UNDER THE THIN LENS APPROXIMATION 

We consider the wave propagation under the gravita- 
tional fields of a lens. We assume that the spacetime 
metric is a Minkowski spacetime plus a small perturba- 
tion due to the existence of a static lensing object. Then 
the metric can be written as 

ds 2 = g^dx^dx" = -(1 + 2U)dt 2 + (1 - 2U)dx 2 , (2) 

where U is a Newtonian potential of the lensing object. 
We consider a propagation of scalar waves <j>, instead of 
gravitational waves, since the wave equation for <f> is the 
same as that for gravitational waves [l^. The scalar 
wave equation, d^^—gg^ d v (f) = 0, with the metric (0) 
is rewritten as 

(A + uj 2 )<j){x) = Aoj 2 U{x)(j){x), (3) 

where we assume that the wave is monochromatic with 
the angular frequency lj. The above equation was 
solved by using the Kirchhoff diffraction integral (see > 
Sec. 4. 7 and 7) under the thin lens approximation. 

We show lensing configuration in Figri] We choose 
a position of the lens as the origin of the polar coor- 
dinate system (r, 9, <ft) . The observer position is r and 

9 = (9 cos (j), 6 sin <p) with 9 <C 1. The incident wave is a 
plane wave propagating in the z-direction. Denoting the 
incident wave as 4>°, we have </>° = e lujrcose ^ 

In this section, we assume the thin lens approximation, 
in which the wave is scattered only on the thin lens plane 
at z = 0, and the lens is characterized by the surface mass 



III. FORMULATION OF NUMERICAL 
CALCULATION 

In this section, we develop a formulation to solve the 
scattering problem of gravitational waves by lensing ob- 
ject which is applicable to the case where the lens po- 
tential is spherically symmetric. A situation we will con- 
sider is that plane wave is entering weak gravitational 
field which is spherically symmetric. When the lensing 
object is spherically symmetric, a scattering problem can 
be reduced to a problem of determining so-called phase 
shift which is used to probe the nature of nuclear physics 
and is also useful for the scattering by BHs where the 
assumption of weak gravitational field breaks down |2l| . 

We have to solve Eq. 10 in order to evaluate quantities 
such as the amplification factor which can be compared 
with the one derived under the thin lens approximation. 
To solve Eq. we choose a center of the lens as the 
origin of coordinate system. Then U becomes a function 
which depends only on the radius coordinate r. Now let 
us expand 4> i n terms of Legendre function 

0(^)=E^^(cos0), (6) 

where £ is a dimensionless variable defined by £ = ru. 
Then equations for <?^(£) are 

where we have used U in stead of U in order to stress 
that U is a function of £. 

For a point mass lens, i.e. [/(£) = — p/(2£), where 
p = 2GMu> and M is a lens mass, analytic solutions of 
Eq. are known as Coulomb wave functions [22^. The 
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solution which is regular at £ = is 

,T(i + l-ip) . t+1 u 



W + i)) 

xF(e + l-ip;2(£ + i);-2i£), 

a e = HTgT(£+l + ip), (8) 

where F is the confluent hypcrgcometric function. The 
solution which is singular at £ = is 



xC/(£ + 1 - ip; 2(* + 1); -2*0 + c.c,(9) 

where U is defined as 
f 



U(a;b;z) 



dte- z H a -\l + t) h ~ a ~\ (10) 



r(a) y 

The asymptotic form of these functions are 

Fl(-P,0 — sin(£+plo g 2£-^-^), 
G*(-p, — cos(£ + p log 2£ - ~* - (11) 

The term p log 2£ in the phase of trigonometric functions 
represents the nature of long range force which is char- 
acteristic of Coulomb force. 

On the contrary, a solution of Eq. (PJ) that a plane 
wave is entering a point mass lens is well known and is 
given by [23| 

4> p = e % p r(l - ip)e^ cos9 F{ip- 1; - cos0)). (12) 

Because the solution Eq. i|12|) is regular at £ = 0, it is 
written as a partial wave sum of regular Coulomb wave 
function, 



t> P = }_^ a e 7 P^(cos0), 



a e = i e (2£+l)e 



-1(7 1 



(13) 



For an extended lensing object, analytic solution of 
Eq. © does not exist usually. However if the lensing 
object exists only in a finite region, then the solution of 
Eq. outside the lens can be written as a summation 
of partial waves which are now a linear combination of 
two independent Coulomb wave functions. By determin- 
ing a coefficient of each Coulomb wave function, we can 
calculate the wave form <p outside the lens. In the aim 
of only determining the wave form far from the lensing 
object, we don't need to know the expression of cj> inside 
the lens. 

Now let us write the solution of Eq. J3J as 



= 4> P + 4>s 



(14) 



Thus <j) s represents the scattered wave which arises due 
to the deviation of lens from a point mass. There may 



be no incoming scattered wave to the lens from infinity, 
so we assume the asymptotic form of </> s as 

00 log 2£ — — 

Mt,0)=J2 —r s^Mcosfl), (15) 



1=0 



2i£ 



where sg are undetermined complex numbers, but not 
arbitrary. In order that tj> in Eq. I|14f> satisfy the wave 
equation J3J), sg must be related to ag as 

Si =a e (e 2iS '-l), (16) 

where 5g are real numbers. In terms of 5g, <fi is written as 



age l5i I cnsi 



+ siufr g ^ ?'® ] Pi(cosfl). (17) 



Thus we can calculate wave form outside the lens object 
by determining the phase shift 5g. Here, Sg are deter- 
mined by matching a solution of Eq. (JJJ with Eq. i|17|) 
at a radius £o being outside the lens. Then, we have 



tan 5/ 



m{tjo)F'g{-p^o) - g'g{tjo)Fg{-p,ti ) 
mo)Gg{-p,t )-gg{i )G' t {-p,t o y 



(18) 



where gg is a solution of Eq. Q which is to be calculated 
numerically. We set the initial condition of Eq. (JJJ at 
£ = is that gg is regular. The range of 5g is from — 7r/2 
to tt/2. 



IV. RESULTS 

We investigate the validity of the thin lens approxima- 
tion for three lens models; the uniform density sphere, 
the singular isothermal sphere and the Hcrnquist model. 



A. Uniform density sphere 

We first present the results for uniform density sphere 
which has the simplest structure next to the point 
mass. The gravitational potential for the uniform density 
sphere is given by 



U(r) 



GM 
' 2R 
GM 



(r < R) 
(r > R), 



(19) 



where r is the distance from the center of the sphere, R 
is the radius of the sphere and M is the lens mass. Here 
we consider the case that R is larger than the Einstein 
radius t\e, in which case the effect of the size of the lens 
object is expected to be important. For R < tei the 
result is almost the same as for point lens mass and it 
is known that amplitude of the amplification factor for 
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the point mass lens coincides with that in the thin lens 
approximation 



The amplification factor in the thin lens approximation 
Eq. © for the case of uniform density sphere is written 
as 



Fthin(r,0) = -i—ei 



dssJ (ujs6y wi ^-^ s) \ 



FIG. 2: h as a function of I for the uniform density sphere. where Jq ig the _ th Begsel function and y is given by 



(20) 



4Mlog l + yjl- (s/R) 2 
AM\og(s/R) 



-iM(4-( S /i?) 2 ) ^\-{s/R) 2 



(s<R) 
(s > R). 



(21) 



r 



Next let us calculate F numerically by using the 
method developed in the previous section. We first show 
5g as a function of I in fig. [21 We see that 8i decreases as I 
increases. In particular, above around I = 90, Si rapidly 
approaches zero. This is because the t -th partial wave 
can be interpreted as an incident particle with impact 
parameter i/ui. Partial waves of I £ ujR pass through 
the lens potential outside the lens object and the grav- 
itational effects on these partial waves are the same as 
point mass, which implies Se becomes zero. 

FigElshows the error of the amplification factor for the 
uniform density sphere as a function of lens parameters 
ujR, wr, GMui and 9 respectively. The vertical axis is 
the error A defined as, A 2 = (|F| - \F thln \) 2 /\F thl . 



The normalized Einstein radius is uite — 
200(GAiw/l) 1 / 2 (wr/10 4 ) 1 / 2 . The radius of sphere loR{= 
100) is comparable to the Einstein radius. 

The top left and right panels show A as a function of 
uiR and uir respectively. We find that averaged in R or 
r over a period of oscillation, A is proportional to R/r. 
Also we find that A remains smaller than R/r which is 
much smaller than unity. (The dashed lines denote the 
ratio of the radius to distance R/r.) This suggests that 
geometrical thickness of lens R/r is a suitable measure 
of the validity of thin lens approximation in wave optics. 

The bottom left panel shows dependence of A on 
GMui. We find that A has a peak around GMui = 10 and 
decreases as GMui increases for GMui > 10. This result 
shows that the thin lens approximation is valid for the 



wavelength where we can use geometric optics instead of 
wave optics. The reason why thin lens approximation be- 
comes valid in the geometric optics is that the deflection 
angle can be evaluated by using the thin lens potential 
if the distance between lens and observer is much larger 
than the lens size (see Ref . [T^| . p. 124). For this panel too, 
A remains smaller than R/r in all range of frequencies 
of calculation. 

The bottom right panel shows dependence of A on 9. 
We find that A becomes maximum at close to 9 = and 
decreases as 9 becomes larger. This is because for large 
9 such that r9 > R (9 > 0.01 for a case of fig. El, the size 
effect of lens becomes negligible and the lensed waveform 
becomes the same as for the point mass lens for which it 
is known that thin lens approximation is valid. Except 
for a small region of 9 around where A takes maximum 
value, A is smaller than R/r. 

To summarize, for the uniform density sphere, the rel- 
ative error of the amplification factor in the thin lens 
approximation is suppressed within the ratio of its ra- 
dius to distance R/r, which is much smaller than unity 
in a realistic astrophysical situations and is the largest 
for A ~ GM (A is the wavelength). 



B. Singular Isothermal Sphere (SIS) model 

We also calculated A for singular isothermal sphere 
which is a model of galaxies, dark matter haloes and 
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FIG. 3: Error of the amplification factor for the uniform density sphere. The horizontal axis is the radius uiR (top left), the 
distance lot (top right), the mass 2GMuj(= p) (bottom left), and the angle 9 (bottom right). The dashed lines denote the ratio 
of radius to distance R/r. 



star clusters. The density profile of SIS model is p{r) = 
v 2 /{2irr 2 ) where v is the velocity dispersion. For nu- 
merical calculation, we have to introduce cutoff radius r c 
because gravitational potential does not approaches to 
tx r 1 far from the lens without cutoff. Thus we here 
assume that the density vanishes for r > r c . Newton 
potential becomes 



where M is the mass inside the cutoff radius r c : M 
2v 2 r r . 



U(r) 



<2M(ln£-l) (r<r c ) 
(r>r c ), 



(22) Fthin is given by Eq.(j2DJ) an d ip(s) is given by 



^s) = I -^f K V^f^+ ^arctan^f^ + AGM In (l + ^l- {s/r s ) 2 ) , (s < r c ) (23) 
{8GMHs/r c ) V (s>r c ) 



Fig. 2| presents dependence of A on parameters, wr c , u>r, ujGM and 0. We find that the qualitative behavior 
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FIG. 4: Same as Fig0 but for the SIS model. The dashed lines denote the ratio of cutoff radius to distance r c /r 



is the same as the case of uniform density sphere. We see 
from the top left panel and the top right panel that thin 
lens approximation becomes worse as cutoff approaches 
to observer's distance r from the center of SIS. 



where r s is a scale length and p s is a characteristic den- 
sity. For numerical calculation, we introduce cutoff ra- 
dius r c for the same reason as SIS. For r < r c the density 
is given in Ea. (|24l) . while for r > r c the density vanishes. 
Then, Newton potential becomes 



C. Hernquist model 

We also calculated A for Hernquist model |l9j, which 
fits well the luminosity distribution of many elliptical 
galaxies and bulges. Its density profile is given by 



p(r) = 



(r/r s )(l + r/r s 



(24) 



GM 

U(r) — < r + r a 

GM 



( r rr k ) 2 + 7: G ' M ( r < r c) 

(r > r c ), 



(25) 



where M is the mass inside the cutoff radius r c : M 
2TT Ps r 3 s [r c /(r s +r c )} 2 . 

Fthin is given by Ea. (|2U[) . and 4>(s) is given by 



— AM Ts+2 n src arccosh— 



-AM r ° + ^ sTc arccosh- 



4M In s + A^My/r 2 - s 2 + 



AM Ins + ApiMy/r 2 



4M Ins. 



4Mr f arctanh 

r c V r 2 —s 2 

r c V S ' 2 ~ V s 



2)( r 2 



s 2 ~\-r c r s 



{s 2 —r 2 )(r 2 —s 2 ) 
s 2 -\-r c r s 



(s < r s ) 

(r s < s < Tcj 
(r > r c ), 



(26) 
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Fig. shows A for Hernquist model. We see that the 
behavior of A is almost the same as SIS. 



V. SUMMARY 

In this paper we discussed the validity of the thin lens 
approximation in the framework of wave optics. 

In sec III, we developed a formalism to solve the wave 
equation for the spherically symmetric potential. In this 
case, the partial differential equation can be reduced to 
a set of ordinary differential equations. The method we 
used is to determine the so-called phase-shift which repre- 
sents the difference of scattered waves between the point 
mass lens and an arbitrary spherically symmetric poten- 
tial. This formalism is only applicable to the spherically 
symmetric lens whose size is finite. For lens models such 
as the SIS profile which extends to infinity, we have to 



introduce cutoff to make the total mass finite. 

We also solved the wave equation numerically for the 
spherically symmetric potential. By numerical calcula- 
tions, we found that the error of the thin lens approxima- 
tion for the simple lens models is the same as or smaller 
than the geometric thickness of the lens, s/r, where s 
is the size of the lens and r is the distance between the 
lens and the observer. The error is the largest for the 
wavelength comparable to the Schwarzschild radius of 
the lens. 
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